Emerging attractors and the transition from dissipative to conservative dynamics 
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The topological structure of basin boundaries plays a fundamental role in the sensitivity to the 
final state in chaotic dynamical systems. Herewith we present a study on the dynamics of dissipative 
systems close to the Hamiltonian limit, emphasising the increasing number of periodic attractors, 
and on the structural changes in their basin boundaries as the dissipation approaches zero. We show 
numerically that a power law with nontrivial exponent describes the growth of the total number of 
periodic attractors as the damping is decreased. We also establish that for small scales the dynamics 
is governed by effective dynamical invariants, whose measure depends not only on the region of the 
phase space, but also on the scale under consideration. Therefore, our results show that the concept 
of effective invariants is also relevant for dissipative systems. 

PACS numbers: 05.45.Ac 61.43. Hv 
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I. INTRODUCTION 

Many dynamical processes have been shown to possess 
coexisting metastable states, and their dynamics can be 
highly sensitive to the initial conditions. Some examples 
include neural behaviour [l], rain events 0], earth- 
quakes' dynamics [J, effective fractal dimension of 
energy levels @, among others. Their correct interpre- 
tation requires an understanding of dynamical systems 
with very low dissipation, lying in between the strongly 
dissipative limit and the conservative one. In spite of 
the importance of this problem, there are few system- 
atic studies of the low-dissipation limit and the transi- 
tion from dissipative to conservative dynamics. This is 
the problem we address in this paper. 

In Hamiltonian systems, chaotic regions typically coex- 
ist with regions of regular motion around the marginally 
stable periodic orbits, also known as Kolmogorov-Arnold- 
Moser (KAM) invariant tori. Chaotic trajectories have 
intermittent behaviour and spend long times sporadically 
near the border of regular islands. This stickiness to the 
KAM tori makes their dynamics fundamentally differ- 
ent from that of hyperbolic chaotic systems, since even 
small islands can exert great influence on the global dy- 
namics @, H, 0- On the other hand, strongly dissipa- 
tive systems are characterised as having only one or few 
attractors to which all the initial conditions eventually 
converge. These two regimes could not be more different 
from each other; however, if the dissipation decreases to 
zero, the phase space structure of the dissipative systems 
must evolve in such a way as to approach the complex 
hierarchical organisation found in Hamiltonian systems. 
How this transformation occurs is far from obvious. 

Despite the existence of extensive investigations on 
strongly dissipative and on Hamiltonian systems, very 
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little is known about properties of systems close to 
the border between dissipative and conservative dy- 
namics, where many coexisting periodic attractors are 
present flfj| |. In this paper we investigate the growth of 
the number of periodic attractors and how the topology 
of the phase space evolves when the dissipation is re- 
duced. We obtain numerically a power law describing 
the growth of the number of periodic attractors as the 
damping approaches zero. When the system is close to 
the Hamiltonian limit, we find that high-period periodic 
attractors become increasingly important in the system's 
dynamics, in contrast to the case with high dissipation. 
These emerging attractors proliferate ever more rapidly 
as the dissipation decreases. Moreover, we argue that 
although this number can be extremely high for small 
dissipation, we expect it to be finite (for non-zero dissipa- 
tion), as conjectured recently [TJ The power law of 
the number of attractors with dissipation, found in this 
work, corroborates this conjecture. Another important 
finding in this paper is that for low levels of damping 
the dynamics is characterised by effective invariants, a 
concept previously used in the context of non-hyperbolic 
dynamics of Hamiltonian systems Q ■ We study specif- 
ically the effective fractal dimension, which we show to 
depend on the scale and also on the region of the phase 
space. Therefore we show that effective dynamics is rele- 
vant also for dissipative systems, and not only for Hamil- 
tonian ones. 

This paper is organised as follows. We start, with 
Sec. HH where we introduce the dynamical systems used 
in our analysis and revisit the problem of the number of 
coexistent periodic attractors. It is followed, Sec. IIIII by 
a discussion regarding the topology of the phase space in 
the low-dissipation limit. In Sec. IIVI we introduce effec- 
tive invariants and define the effective fractal dimension 
for weakly dissipative systems. Finally, Sec. IVl brings our 
conclusions. In the appendix, we bring up some formal 
definitions. In particular, we recall a definition of periodic 
attractors, and discuss its implications for both experi- 



mental and numerical realistic investigations, where one 
typically is only able to make use of limited precision. 
We also recall a definition of basin of attraction, which is 
coherent with the definition of periodic attractors used 
here. 



II. NUMBER OF PERIODIC ATTRACTORS 

A. Single Rotor Map 

A paradigmatic system that allows the investigation 
of the transition from the dissipative to the conservative 
limit is the single rotor map, obtained from the mechan- 
ical pendulum kicked periodically at times nT, n G N + , 
with force fo [U[: 

x n +i = x n + y n (mod2Tr) 

y n+ i = (1 - v)y n + f sin(x„ + y n ), 

where x G [0, 27r] corresponds to the phase, y G R to the 
angular velocity, and v G [0, 1] is the damping parameter. 
For v = 0, the well known area preserving standard map 
analysed by Chirikov [l4| is recovered. For v = 1, we 
obtain the circle map with zero rotation number. 

Previous work has considered the question of the num- 
ber of attractors for this system They have derived 
an analytical expression for the number of period- 1 pri- 
mary attractors, that have not undergone period dou- 
bling bifurcation: 

^=2/(^) + l, (2) 

where /(■) denotes the integer part of the expression in 
brackets. As pointed out in Ref. [l(|, numerical detec- 
tion of attractors with high period is a difficult task. The 
reason is the extremely small size of their basins of at- 
traction, what requires computation using very fine grids. 
Adding to this, they have short lifetime in the space of 
parameters. Therefore it is necessary to follow a great 
number of initial conditions in order to be able to detect 
them. 

Here we overcome this problem. Instead of choosing 
an a priori fixed number of initial conditions in a preset 
grid, as previously done in Ref. [h| and elsewhere, we fix 
the non linearity fo = 4.0 for a given damping v and then 
iterate an ensemble of randomly chosen initial conditions 
with uniform distribution on a bounded region of the 
phase space. Therefore, the chance for an attractor to 
be found does not depend on our particular choice of the 
grid. For the dissipative case the dynamics takes place 
in the cylinder T = [0, 2-7r] x R and the area where the 
initial conditions are chosen depends on the damping. 
Since v > 0, from the second equation in the map (fTJ), 
one gets |y fc+1 | < (1 - v)\y k \ + f , so if \y k \ > f Jv => 
\Vk+i\ < \Vk\- Hence the attractors are found in A C T, 
with A = [0, 2ir] x [-y ma x,ymax], where y max = f /u. 
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FIG. 1: Total number of periodic attractors detected for a 
specific value of the dissipation v, as a function of the num- 
ber of initial conditions randomly chosen in A and evolved 
according to the map. The parameters for the map {1} were 
fo = 4.0, and v = 0.0035. 
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FIG. 2: (Colour online) Total number of periodic attractors 
which were detected numerically for the single rotor map 
(black circles), and for the Henon map (red diamonds) for 
different values of v. The parameters were fo = 4.0, for the 
single rotor map, and A = 1.075 for the Henon map. It is also 
shown the power laws fits, the black and the red solid lines 
for the single rotor map and for the Henon map, respectively; 
compare to the number of period 1 attractors expected from 
Eq. ([2]) (black dashed line) for the single rotor map. The 
divergence from this curve (dashed black line) to the other 
(black solid line) is clear, showing the important contribu- 
tion of high period periodic attractors in the regime of small 
damping for the single rotor map. 



In order to detect the attractors, even with very small 
basins of attraction, we keep iterating new, randomly 
chosen initial conditions (xo,yo) G A until their trajec- 
tories converge to some periodic attractor. We use up 
to 2 x 10 5 iterations; independently of the damping, this 
provides a convergence of more than 99.5% of the initial 
conditions. We have ignored the remaining trajectories 
(which may include chaotic orbits), regarding them as 
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not being statistically representative. It is known that 
for this and similar systems, almost all attractors consist 
of periodic orbits. We keep iterating new initial condi- 
tions until we fail to find new attractors. This was deter- 
mined by the criterion that whenever the total number 
of detected attractors did not change for the last 10 6 
initial conditions, we assume that all the detectable at- 
tractors have been found. This is illustrated in Fig. [1] for 
v = 0.0035. We repeat this procedure for different values 
of v 1 in order to find how the total number of periodic 
attractors changes with damping. Notice that the area of 
the phase space where the trajectories are trapped grows 
as we decrease v; for v = 0.001 for example, we need 
22 x 10 6 initial conditions to find most attractors in the 
system. 

The question regarding the finiteness of the number of 
attractors and their density in phase space has been con- 
sidered before and it is widely regarded as one of the most 
important open problems still to be answered in Dynam- 
ical Systems Theory [Tl| - Among many conjectures, it 
has been initiall y p roposed that the number of attrac- 
tors was infinite [l5j . Afterwards it has been proved that 
this should only hold for a zero measure set of parameter 
values, though dense in some interval Recently it 

has also been conjectured by Palis that the total number 
of attractors should be finite [l2j . 

To investigate this issue using our system (JTJ) as a test- 
ing ground, we calculate the number of attractors as de- 
scribed above, and see how this depends on the dissipa- 
tion v as it is decreased and approaches the conservative 
limit, i.e. v — ► s > 0. Figure [2] shows a log-log plot of our 
data, yielding a power law describing the dependence of 
the total number of periodic attractors on v. 
Numerical fitting gives us the following law: 



N- 



TP„ 



1.26' 



-1.21 



(3) 



where Ntp r is the total number of detected periodic at- 
tractors for the map (Q]). 

Now we want to compare our power law for the growth 
of the total number of periodic attractors with the for- 
mula © [l(| for the number of period 1 attractors. We 
observe in Fig. [5] that when the damping is decreased 
the contribution of higher period periodic attractors be- 
comes important for the total number of detected attrac- 
tors. The main reason is that the lifetime of a stable, 
periodic orbit in parameter space increases as damping 
is reduced. By lifetime of an orbit we mean the range 
in parameter space for which the orbit exists, without 
undergoing any bifurcation. Since the higher the period 
of an attractor, the smaller its lifetime interval in pa- 
rameter space (Toj . the extension of lifetime for small 
damping promotes the overlap of different periodic or- 
bits in the space of parameters, which was not possible 
before. Equation © is only valid within the parameter 
intervals where none of the period 1 attractors has un- 
dergone period doubling [Io| . Since some of the period 
1 attractors may have gone through such process for the 
parameters in Fig. [21 the difference between the expected 



number of period 1 attractors and the total number of 
attractors can be even larger than that inferred from the 
above argument. 



B. The Henon Map 

We further investigated the growth of the number of 
periodic attractors for a map of a different class than the 
previous one. We chose the Henon map in the form [13]: 



x n +\ = A- x n - (1- v)y n 

i/n+l — ^n; 



(4) 



where A represents the bifurcation or non linearity con- 
trol parameter. The parameter v £ [0, 1] represents the 
damping parameter. When v = 1, the two equations in 
Eq. [5] are no longer coupled, and we obtain the quadratic 
map. For the other limit, v = 0, we have a conserva- 
tive map, hence the determinant of its Jacobian matrix 
is equal to 1. Contrary to the case for the map |T]), the 
dynamics for the Henon map Q is not a priori contained 
within a region of the phase space. For a range of param- 
eters, the initial conditions can either be trapped by some 
of the coexistent periodic attractors, or be scattered. In 
fact, most of the initial conditions diverge to infinity. 

We repeated the procedure applied to the previous 
map in order to obtain the number of periodic attrac- 
tors for the Henon map. We fixed A = 1.075 for 
the map ([J), and we iterated a set of initial conditions 
in [—5, 5] x [—5, 5] verifying whether they converged to 
some periodic attractor or diverged. Because the num- 
ber of initial conditions converging to periodic attractors 
is much smaller than that diverging, in order to obtain 
about 20 x 10 6 initial conditions converging to periodic 
motion, as many as 20 x 10 7 initial conditions were nec- 
essary depending on the damping. 

Figure [2] shows a log-log plot of our data, yielding a 
power law describing the dependence of the total number 
of periodic attractors on v. For the Henon map, the 
numerical fitting produces: 



N- 



tp h 



0.025 ■ 



-1.61 



(5) 



where Ntp h is the total number of detected periodic at- 
tractors for the map (g]). 

As previously observed 17[, the Henon map belongs 
to a class of dynamical systems whose the conservative 
element has few primary island, which are surrounded by 
secondary ones. Therefore, although it also observed the 
coexistence of periodic attractors for a range of param- 
eter, it is expected this number to be smaller than that 
for the single rotor's family. 



C. On the General Behaviour 

Even though we cannot claim the above result to be 
a mathematical proof, our results seem to uphold Palis' 
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conjectures [HI, [11], i- e -> f° r any arbitrarily small e > 0, 
the limit of Eqs. ©, [©], as lim„_ e + N TPr[h] predicts 
a finite number of attractors for the system, though this 
number can be very high for small e (Ntp r , and Ntp h 
for the single rotor map, and for the Henon map, respec- 
tively). 

Although some attractors may have their basins of at- 
traction so tiny that even by iterating a huge number 
of initial conditions is not enough to numerically detect 
them. Figure [1] shows that the number of attractors 
which was detected is asymptotic to a constant (we ob- 
tain a similar figure for the Henon map). Therefore, our 
results indicate that the number of attractors increases as 
a power law as the dissipation decreases, being finite for 
any given v > 0, though arbitrarily large for v approach- 
ing zero. Even though the specific value of the exponent 
in Eq. ((3|) and Eq. (J5]) may be only valid for the map (fT]), 
and the map (0]), respectively, the growth of the number 
of periodic attractors is a very generic property of weakly 
dissipative systems [TtJ which should not depend on the 
details of the system, since in the low-dissipativc limit 
the system's phase-space approaches that of Hamiltonian 
systems, which have universal properties. We therefore 
expect that the same behaviour to hold for other systems. 



The natural question which arises is the following: 
what happens to the topology of the phase space as new 
attractors continually appear with decreasing damping? 



III. TOPOLOGY OF THE PHASE SPACE 

We start this section by presenting a qualitative picture 
of the transformations occurring in phase space struc- 
ture when the damping is decreased. For simplicity, we 
restrict ourselves to the map (U}. Although one cannot 
draw definitive conclusions from this qualitative analysis, 
it motivates us to some further quantitative investigation, 
which will be tackled at the end of this section and in the 
next one. In order to get some clue about which topo- 
logical changes occur when the damping is decreased, we 
begin by looking at the phase space for the system with 
constant forcing and high damping. In this case the sys- 
tem has only few periodic attractors. We decrease then 
i>, find the periodic attractors for the new parameter and 
plot their basins of attraction. Since we are now only 
interested in the general picture, the number of initial 
conditions is fixed at 2 x 10 6 . 



If we follow the distributions of filaments of the largest 
basin of attraction in the phase space, we notice that 
their structure becomes more heterogeneous for different 
regions when the damping is decreased. This is illus- 
trated in Fig. [3l where each colour (colour online) rep- 
resents a set of initial conditions which converge to one 
of the attractors, i.e., different colours represent different 
basins of attraction. In Fig. [3] (a), for which v = 0.32, 
we notice that the largest basin of attraction (in black) 
is spread out over the phase space in an intricate struc- 
ture. Moreover, if we compare distinct regions of the 
phase space, there is a qualitative difference in the "den- 
sity" of filaments of the largest basin of attraction. In 
particular, the farther the region is from the main at- 
tractor, the thinner the filaments become. In Fig. [3] (b), 
for v = 0.3, we detect five attractors and the topology 
of the phase space, i.e., the spatial distribution of invari- 
ant sets, seems to become more complex. Going even 
further, for v = 0.2415 in Fig. [3] (c), we have seven at- 
tractors in the same region and the difference of density 
of those filaments is clearer. Figure [3] (d) is the blow up 
of the region x E [4,5], G [-15,-10] for v = 0.2415. 
The rescaling of the phase space shows the structure of 
the basins which seems to form a self similar pattern, as 
illustrated in Fig. |U 

As we decrease the damping even further, this pro- 
cess seems to carry on and the heterogeneity of the 
phase space keeps increasing. Furthermore, the differ- 
ence in density of filaments of distinct regions of the 
phase space is seen by looking, for instance, at the largest 



basin of attraction for v = 0.02, and fo = 4.0. For 
these parameters we expect from Eq. to find about 
140 periodic attractors. The largest basin of attraction 
is responsible for about 36% of the initial conditions 
(xo,Vo) £ [0, 27r] x [— 7r,7r] [Icj . This is illustrated in 
Fig. [5l which has been done by the iteration of an en- 
semble of 10 6 initial conditions. 

The structures we just presented are nothing but man- 
ifestations of some of the underlying invariant sets. That 
is, the dynamics has multiple periodic attractors, and the 
closure of the set of initial conditions that approaches 
a give attractor is its basin of attraction. The bound- 
ary which separates different basins of attraction is the 
basin boundary. This basin boundary is the stable man- 
ifold of an invariant set. If we define the boundary B in 
some region of the phase space embedded in the bound- 
ary as S n B, where S denotes some open set in phase 
space which contains part of the boundary B, we ex- 
pect for a given v the fractal dimension of the basin 
boundary, dim(S n B), to be constant at different re- 
gions of the phase space, R C I\ as it has been conjec- 
tured in Ref. [l8j]. Nevertheless, for physically relevant 
scales, the heterogeneity in the distribution of the basins 
over the phase space, and hence the spatially heteroge- 
neous distribution of invariant manifolds, suggests that 
course-grained measurements of quantities such as the 
fractal dimension should lead to different results depend- 
ing on the region of the phase space and on the scale 
under consideration. The reason is that realistic mea- 
surements never take the limit of infinitely small scales, 
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FIG. 3: (Colour online) Basins of attraction for fo — 4.0 at different values of v. Each colour represents a set of initial conditions 
which converge to one of the attractors. In (a) v — 0.32, in (b) v = 0.3, in (c) v = 0.2415, and (d) is the blow up of a region 
x G [4, 5], y £ [—15, —10] for v = 0.2415. We used 2 x 10 6 initial conditions for each one of the figures. 



but must have a finite lower scale. Previous works in 
the context of Hamiltonian systems show that when the 
phase space has a heterogeneous structure such as in this 




FIG. 4: (Colour online) Blow up of the region x £ 
[4.662,4.679],?/ <E [-13.955,-13.94] for v = 0.2415 in Fig. ES 
(d), showing the fractal Cantor structure. 



case, one needs to go to exceedingly small scales to ap- 
proach the mathematical value of the fractal dimension, 
which is unique and independent of the portion of the 
phase space used to calculate it. For physically realistic 
scales, an approximation to the fractal dimension (and 
similar quantities) is more appropriate to describe the 
system's dynamics; this approximation, called the effec- 
tive fractal dimension, depends on the position as well as 
on the scale under consideration. This effective fractal 
dimension will be presented in more detail in the next 
section. 

A finite-scale approximation of the fractal dimension 
of the basin boundary can be estimated using the uncer- 
tainty method [l9j . It consists in iterating an ensemble of 
initial conditions (xo,yo) and checking to which attrac- 
tor they converge. Then we add a small perturbation 
e to every initial condition, say (xo,yo + £ )' an< ^ check 
whether they change from one basin to another. We 
count the fraction of initial conditions that have changed 
their asymptotic state after being perturbed in this man- 
ner. This fraction, f(e) of e-uncertain points is the frac- 
tion of initial conditions that change basins under an e- 
size perturbation, and it scales as /(e) ~ e a , where a 



6 



Ill III 



FIG. 5: The largest basin of attraction for v = 0.02 and /o = 
4.0. The black points represent initial conditions which con- 
verge to the attractor formed from the main island. The initial 
conditions were chosen such that, (xo,yo) £ [0, 2ir] x [— n, n]. 
We used the same parameters as those used in Ref [lOj ] . 




[1.0,2.5] x [198.5,200.0], and they were chosen with the 
same Lebesgue measure. We chose three different regions 
of the phase space, located at different distances from the 
main attractor. We remark that, although such regions 
may be out of A C T defined in section [HI it does not 
follow that they do not contain boundaries. This is be- 
cause regions A C T, defined for each value of v, are 
portions of the phase space where the attractors are to 
be found; their basins, on the other hand, are expected to 
be extended throughout the phase space. Since the den- 
sity of periodic attractors in phase space decreases when 
we look at regions farther away from the main attrac- 
tor, we expect the stable manifold to be denser around 
the central part of the phase space. Furthermore, given 
that we expect the dissipative regime to approach the 
complex structures of the Hamiltonian dynamics at the 
zero dissipation limit, it is sensible to expect a distribu- 
tion of manifolds which are similar for conservative and 
for low dissipative dynamics. For Hamiltonian dynamics, 
the density of manifolds decreases as one looks at regions 
farther from the main island Q- 

Looking at Fig. [IJJ it is clear that for e > 10~ 13 , the 
estimated fractal dimension is larger for the region R3, 
where we expect the stable manifold to be less dense. 
Only when we go down to e < 10~ 13 does the slope for 
the region R3 converges to the same value as the one 
of the region R±; for small enough e, the slope eventu- 
ally converges to the true fractal dimension, which has 
a unique value [Isl ]. We have computed /(e) for even 
smaller values of u, and we have observed the same be- 
haviour. However, for values as small as v = 0.02, even 
going down to e = 10 -35 is not enough to obtain a con- 
vergence to a unique value of the estimated dimension in 
different regions. This shows that an effective dimension 
is indeed the relevant physical quantity to be considered 
also for weakly dissipative systems. 



FIG. 6: Fraction of uncertain initial conditions /(e) as a 
function of the perturbation size e for two different regions 
of the phase space. The regions, Ri C T, were defined as 
Rx = (x,y) G [0,1.5] x [-1.0,-2.5], and R 3 = (x,y) € 
[1.0,2.5] x [198.5,200.0]. The utilised parameters were f = 
4.0, and v = 0.08. 



is related to the box counting dimension d of the basin 
boundary by a = D — d, where D is the dimension of the 
phase space (For the map ([T]) we have D = 2). 

In Fig. [5] we show the fraction of "uncertain" initial 
conditions as function of the size of perturbation e [2(| • 
The slope of the curves are related to the fractal dimen- 
sion of the basin boundary. 

The exponents a = Aln/(e)/Alne can be computed 
over some decades. However, we notice in Fig. [6] that 
the slope assumes different values for distinct regions of 
the phase space. Such regions, Ri C I\ are defined as 
J?i = (x,y) G [0,1.5] x [-1.0,-2.5], and R 3 = (x,y) G 



IV. EFFECTIVE FRACTAL DIMENSION 

From the theoretical results regarding the fractal di- 
mension of the basin boundary (l8| , we expect to have the 
same fractal dimension for different regions of the phase 
space. However, the above discussion and results suggest 
that, for realistic scales, the behaviour for weakly dissipa- 
tive systems can be quite different from the asymptotic 
one. Therefore it is sensible to describe the system in 
terms of effective invariants, a concept which has been 
used mostly in the context of non-hyperbolic Hamilto- 
nian dynamical systems. 0]. As an example of effective 
invariants we use the definition of effective fractal dimen- 
sion for a D-dimcnsional system, which depends on the 
region R C T, and on the scale under consideration [22| . 
It is defined as 



D eff (R,s)= D- 



dln/(e') 



d\ne' 



(6) 
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FIG. 7: Fraction of uncertain initial conditions /(e) scaling 
with the perturbation e for three different regions of the phase 
space. The regions, Ri C T, were defined as Ri = (x,y) £ 
[0,1.5] x [-1.0,-2.5], R 2 = (x,y) £ [1.0,2.5] x [18.5,20.0], 
and R 3 = (x,y) £ [1.0,2.5] x [198.5,200.0]. The utilised 
parameters were /o = 4.0, and v = 0.3. 




FIG. 8: Fraction of uncertain initial conditions /(e) scaling 
with the perturbation e for three different regions of the phase 
space The regions are the same as in Fig[7] The utilised pa- 
rameters were /o = 4.0, and v = 0.07. 



FIG. 9: Fraction of uncertain initial conditions /(e) scaling 
with the perturbation e for three different regions of the phase 
space. The regions are the same as in Fig[7] The utilised 
parameters were /o = 4.0, and v — 0.02. 



apparently 5 periodic attractors. We notice that for R3, 
where the stable manifold is less dense, starting from 
e = 1CP 1 , the value of the effective fractal dimension is 
larger, but after few decades it converges to nearly the 
same as for R\, i.e., D e ft = 1.89, where the stable man- 
ifold is denser. The same phenomenon is observed for 
i?2, though the convergence is faster than for R3. When 
we decrease v, hence increasing the number of attrac- 
tors, even for e = 10~ 12 , which is almost in the limit 
of normal computation and far beyond realistic measure- 
ment capacity, we do observe different values of effective 
fractal dimension for different regions, as it is shown in 
Fig. [5]and Fig. [HI for v = 0.07 and v = 0.02, respectively. 
We also notice that for a given region the exponent a 
becomes smaller as the damping is decreased. This indi- 
cates that for such weakly dissipative systems, the frac- 
tal dimension is very close to the dimension of the phase 
space. This is, of course, in accordance with the fact that 
for Hamiltonian systems [y = 0), the fractal dimension 
assumes its maximum value, i.e., D. 



CONCLUSION 



where /(e') is the total uncertain phase space volume 
estimated as e' D N(e'), and N(e') is the number of hy- 
percubes of size e' necessary to cover S (~l B. 

We have computed the effective fractal dimension of 
the basin boundary for finite scales, for different val- 
ues of is, and for distinct regions in the phase space. 
The regions, Ri C T, were chosen to be R\ = (x, y) £ 
[0, 1.5] x [-1.0, -2.5], R 2 = (x,y) £ [1.0, 2.5] x [18.5, 20.0], 
and # 3 = (x,y) £ [1.0,2.5] x [198.5,200.0]. 

We have extracted the effective fractal dimension from 
Fig. [7] for v — 0.3 and fo = 4.0, where the system has 



We have shown that the dynamics of weakly dissipa- 
tive dynamical systems can be quite different from cither 
the strongly dissipative systems or the Hamiltonian ones. 
When the damping is decreased, the number of periodic 
attractors increases. We have shown that the growth of 
the total number of attractors is described by a power 
law. In general, the interval in the space of parameters 
in which an attractor exists is smaller for high period 
periodic attractors. For small damping, attractors with 
different periods coexist in phase space, and the contribu- 
tion of high period periodic attractors to the dynamics 
becomes important. Although the total number of at- 
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tractors in the system is very high for low damping, our 
results strongly suggests that it remains finite for v > 0, 
supporting Palis' conjecture. The formation of new at- 
tractors as the clamping decreases changes the topology 
of the phase space considerably. For low damping, the 
dynamics is better characterised by effective dynamical 
invariants, in a similar way to the case of non-hyperbolic 
Hamiltonian systems. These effective dynamical invari- 
ant sets depend on the scale and they differ from one 
region to the other in the phase space. In particular, 
the effective fractal dimension of the basin boundary is 
larger for regions where the stable manifold is less dense. 
For a given region, the smaller the damping, the closer 
the effective fractal dimension is to the dimension of the 
phase space. Although we have used a specific example 
to illustrate our ideas, our results are generic since we 
have deal only with general characteristics of dynamical 
systems close to Hamiltonian case [l2[ ■ 
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Appendix 

On Periodic Attractors 

Some of the most important characteristics of the dy- 
namics are related to invariant subsets in the phase space, 
which attract their neighbouring points. These sets are 
called attracting sets, or attractors, depending on the con- 
text. There are various definitions for attractor, the main 
difference is related to which points in the neighbourhood 
must approach the set. We use the definition given for a 
diffeomorphism f : M. N — > M. N , that is / and its inverse 
arc diffcrcntiablc and the partial derivatives are continu- 
ous. The main source of our definition is Ref. [23|. 

The dynamics takes place in a phase space T C M. N , 
where there is the notion of distance disi(xi, X2) be- 
tween points Xi and X2, both belonging to T. For ex- 
ample, G?is£(xi, X2) can be the Euclidean distance [24| . 
We say that a bounded and limited region of the phase 
space, i.e. Ua C T is called a trapping region for / if 
/(Ua) C int(UA), where int(UA) stands for the inte- 
rior of Ua- Another important concept for the defini- 
tion of attractors is the notion of e- chain. An e- chain 
of length n from x to y, for a map / is a sequence 
{x = xo,...,x„ = y}, such that for all 1 < j < n, 



we have dist(f(x.j-i),Xj) < e. Now we are ready to 




FIG. 10: Illustration of the dependence of the number of peri- 
odic attractors with the precision. The figure shows attractor 
Xi, and yi, with i = 1, or 2. Hence, the attractors are of 
period two. If on the detection of such attractors, e > £; is 
used, only one attractor of period 2 is detected. On the other 
hand, for e < £, we are able to distinguish them. 



state a formal definition of periodic attractors. If for- 
ward iterating our map / in a trapping region, i.e. if 
we take the set A = {~\ n>0 ,f n (U a) , an< ^ we eventually 
end up in an unique closed sequence of points such as 
these points are an e- chain no matter how small e > 
is, we have precisely what is defined as attractor. In 
some other words, a periodic attractor is the unique se- 
quence of points {xi, X2, . . . , x p , xi} remaining by iter- 
ating a trapping region, such that disi(/(xj_i), Xj) < e, 
for all e > 0. In this case, we have period p. Notice that 
the attractivity of the attractor is ensured by requiring 
the existence of the trapping region. Furthermore, we 
remark that, finding multiple periodic attractors require 
the coexistence of different trapping regions. 

Although the formal definition of attractor requires the 
existence of such sequence of points, the e- chain, for all 
e > 0, and infinitely many iterations, both experimen- 
tally, and numerically, we are constrained by limited pre- 
cision. Hence we are forced to use a mild condition on e. 
Therefore, the number of periodic attractors obtained in 
realistic investigation intrinsically depends on the preci- 
sion that we choose. We give a picturesque illustration 
of it in Fig 1 101 On determining of the number of periodic 
attractors in this paper, we have used e = 10~ 10 . Nev- 
ertheless, it does not minimise the importance of ours 
findings, showing fundamental characteristics of the dy- 
namics within realistic scale. 

We also use a more general definition of basin of at- 
traction. The basin of attraction of some trapping re- 
gion U is the set of positive Lebesgue measure of initial 
conditions, whose orbits eventually enter U , as defined 
in Ref. [25|. Note that this slightly different definition 
avoids problems of coherence for not having an infinitely 
small e. Because any periodic attractors is contained 
in some trapping region, hereafter we often mention the 
basin of attraction of some attractor, and it is under- 
stood the basin of a given trapping region containing the 
attractor under consideration. 
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